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EXECUTIVE  SUMMARY 


A.  OBJECTIVE 

The  objective  of  this  research  effort  was  to  assess  the  accuracy  of  using  the  shock  response 
spectrum  (SRS)  for  defining  equipment  shock  tolerance.  The  focus  of  the  research  was  on 
mission>critical  equipment  installed  in  airbase  hardened  structures,  designed  to  survive  the  airblast 
and  ground  shock  effects  of  a  close-in  detonation  by  a  modem,  nonnuclear  weapon.  The 
motivation  was  the  prospect  of  improving  the  shock  tolerance  of  such  equipment,  by  first 
improving  equipment  shock  tolerance  characterization,  thereby  improving  knowledge  of  the 
factors  controlling  equipment  shock  tolerance. 

B.  BACKGROUND 

A  shock  response  spectrum  is  a  plot  of  the  maximum  response  amplitude  of  a  damped  single- 
degree-of-ffeedom  (SDoF)  oscillator  subjected  to  support  excitation,  as  a  function  of  the  natural 
fi-equency  and  damping  value  of  the  oscillator.  An  SRS  is  a  convenient  means  of  illustrating  the 
peak  response  of  a  collection  of  SDoF  systems  to  a  particular  shock  excitation  Due  to  its 
simplicity,  the  SRS  has  become  widely  employed  as  a  means  of  describing  the  shock  response  of 
structures  and  equipment. 

Since  in-structure  motions  are  typically  described  using  shock  response  spectra  (actually 
tripartite  shock  response  spectra),  it  was  only  natural  to  attempt  to  quantify  equipment  failure  in 
terms  of  shock  spectra  The  advantage  of  expressing  equipment  failure  with  the  same 
representation  used  to  quantify  the  in-structure  shock  environment  is  obvious.  Expressing 
equipment  fragility  in  terms  of  shock  response  spectra  greatly  simplifies  the  design  of  equipment 
shock  isolation. 

The  approach  employed  to  determine  equipment  fragility  is  to  subject  an  item  of  equipment  to 
a  base  motion  of  a  certain  intensity,  frequency,  etc  If  the  equipment  does  not  fail,  the  intensity  of 
the  base  motion  is  increased  until  failure  occurs  The  shock  response  spectrum  of  the  base  motion 
which  just  causes  the  equipment  to  fail  is  called  the  shock  tolerance  spectrum,  or  the  fragility 
spectrum  of  the  equipment,  even  though  the  test  SRS  is  input  waveform-dependent 

By  comparing  the  fragility  spectrum  of  an  item  of  equipment  to  the  in-structure  SRS  at  the 
proposed  equipment  location,  shock  isolation  requirements  are  easily  assessed  If  the  equipment 
fragility  spectrum  exceeds  the  in-structure  SRS  at  all  frequencies,  no  shock  isolation  is  required. 


If  the  in-structure  SRS  exceeds  the  equipment  fragility  spectrum  at  any  frequency,  shock  isolation 
is  required. 

This  approach  to  shock  isolation  assessment  is  oversimplified,  and  ignores  several  sources  of 
error.  The  most  important  source  of  error  is  the  fact  that  an  SRS  does  not  correspond  to  a  unique 
input  time-history.  An  infinite  number  of  different  base  motions  can  generate  a  given  SRS.  These 
different  base  motions  could  vary  greatly  in  duration,  frequency  content,  and  amplitude.  The  main 
assumption  behind  the  SRS  approach  to  equipment  fragility  is  that  equipment  failure  is 
independent  of  the  input  waveform.  It  is  assumed  that  the  frequency  content  of  the  support 
motion  does  not  play  a  role  in  equipment  failure.  All  input  base  motions  are  assumed  to  result  in 
the  same  failure  mode.  The  possibility  of  a  single  item  of  equipment  possessing  multiple  failure 
modes  is  not  considered  In  reality,  equipment  fragility  spectra  are  probably  only  valid  for 
frequencies  close  to  the  frequency  at  which  the  item  of  equipment  was  actually  tested 
Extrapolating  equipment  fragility  based  on  nuclear  testing  to  the  conventional  weapon 
environn.ent  is  now  done  routinely.  The  validity  of  this  extrapolation  has  never  been  verified. 

C  SCOPE 

This  report  quantifies  some  of  the  errors  introduced  by  assessing  the  shock  tolerance  of 
equipment  via  the  classic,  linear.  SDoF  shock  response  spectrum.  Two  simple  mathematical 
equipment  models  were  developed  for  the  investigation:  (1)  a  spring  pendulum,  which  possesses 
two  degrees  of  freedom  and  geometric  nonluiearity,  and  (2)  a  clamped  beam  carrying  a  tip  mass 
having  both  translational  and  rotary  inertia,  which  possesses  an  infinite  number  of  linear, 
orthogonal  modes  These  models  were  used  to  investigate  the  effect  of  factors  such  as: 
waveform  dependency  of  equipment  shock  response,  multi-degree-of-freedom  (NlDoF)  vs  SDoF 
equipment  response,  the  effect  of  multidirectional  loading  or  equipment  response,  and  the  effect 
of  geometric  nonlinearity  on  equipment  response. 

D  CONCLUSIONS 

This  study  uncovered  and  quantified  several  shortcomings  of  the  SRS  approach  to 
characterizing  equipment  shock  tolerance: 

•  It  was  proven  that  SRS-based  fragility  spectra  are  not  unique  Even  for  an  item  of  equipment 
which  can  be  modeled  by  a  simple  SDoF  undamped  oscillator,  each  different  base  excitation 
produces  a  different  fragility  spectrum 

•  The  spring  pendulum  model  was  used  to  illustrate  the  response  of  a  geometrically  nonlinear 
item  of  equipment  to  two-dimensional  base  motion  It  was  shown  that  extiemely  erratic 
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behavior  is  possible  if  the  frequency  of  the  base  motion  approaches  one  of  the  natural 
frequencies  of  the  equipment  model.  The  maximum  response  can  be  several  times  greater 
than  that  of  a  simple  SDoF  oscillator. 

•  The  cantilever  beam  model  was  used  to  show  that  the  maximum  repsonse  of  an  MDoF  system 
will  always  exceed  the  response  of  a  simple  oscillator 

In  conclusion,  this  study  has  proven  the  inadequacy  of  the  linear,  SDoF  shock  response 
spectrum  for  characterizing  the  shock  tolerance  of  equipment  Complex,  nonlinear  mechanical 
equipment  subjected  to  multi-directional  support  motion  cannot  be  adequately  represented  by  an 
SDoF  SRS 

E.  RECOMMENDATIONS 


Based  on  the  results  of  this  study,  a  more  rigorous  approach  for  assessing  equipment  shock 
tolerance  is  required.  The  following  points  should  be  kept  in  mind  while  developing  this  new 
procedure: 

•  When  an  item  of  equipment  is  tested  to  determine  its  shock  tolerance,  the  test  input  waveform 
must  be  representative  of  the  anticipated  threat  input  waveform.  Multi-directional  support 
motion  must  be  reproduced.  Equipment  tests  should  excite  the  same  response  modes,  and 
produce  the  same  failure  modes,  as  the  actual  in-service  base  motion. 

•  Analytical  equipment  models  must  be  detailed  enough  to  reproduce  the  salient  features  of  the 
actual  equipment  response.  The  model  need  not  encompass  the  entire  item  of  equipment,  but 
it  must  adequately  represent  the  critical  components  The  input  motion  to  the  model  must 
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SECTION  I 


INTRODUCTION 

A.  OBJECTIVE 

The  objective  of  this  research  effort  was  to  assess  the  accuracy  of  using  the  shock  response 
spectrum  (SRS)  for  defining  equipment  shock  tolerance  The  focus  of  the  research  was  on 
mission-critical  equipment  installed  in  airbase  hardened  structures,  designed  to  survive  the  airblast 
and  ground  shock  effects  of  a  close-in  detonation  by  a  modem,  nonnuclear  weapon.  The 
motivation  was  the  prospect  of  improving  the  shock  tolerance  of  such  equipment,  by  first 
improving  equipment  shock  tolerance  characterization,  thereby  improving  knowledge  of  the 
factors  controlling  equipment  shock  tolerance. 

B  BACKGROUND 

A  shock  response  spectrum  is  a  plot  of  the  maximum  response  amplitude  of  a  damped  single- 
degree-of-freedom  (SDoF)  oscillator  subjected  to  support  excitation,  as  a  function  of  the  natural 
frequency  and  damping  value  of  the  oscillator.  For  the  SDoF  oscillator  shown  in  Figure  1 ,  the 
response  is  controlled  by  three  parameters:  the  mass,  w;  spring  stiffness,  k,  and  damping 
coefficient,  c  If  the  base  motion  (i.e.,  shock  excitation)  is  described  by  a  displacement  time- 
history,  y{t),  and  the  resulting  displacement  of  the  mass  is  given  by  x(/);  the  relative  displacement 
between  the  mass  and  the  base  may  be  calculated  as,  w(/)  =  x(r) -><(/)  The  peak  SDoF  system 
responses  of  most  interest  in  shock  analysis  are;  the  maximum  absolute  value  of  relative 
displacement  between  the  mass  and  base  (spectral  displacement,  the  maximum  absolute 
value  of  relative  velocity  between  the  mass  and  base  (approximated  by  the  spectral  velocity,  i',.), 
and  the  maximum  absolute  value  of  mass  acceleration  (approximated  by  the  spectral 
acceleration,  S^).  The  three  spectral  quantities  are  related  by  the  following  defining  equation 

— =  (1) 

where 

03„  =  J—  (rad/sec)  (2) 

V  m 

is  the  natural  circular  frequency  of  the  SDoF  system  This  relationship  allows  spectral 
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igure  1  Smgle-uegree- 


displacement,  velocity,  and  acceleration  to  be  plotted  on  a  single  graph,  as  a  function  of  system 
natural  circular  frequency,  O),;  or  system  natural  frequency,  /„;  where 

(Hz)  (3) 

Such  a  plot  is  referred  to  as  a  tripartite  shock  response  spectrum.  Typical  tripartite  shock 
response  spectra  are  shown  in  Figure  2.  It  is  common  for  a  tripartite  shock  response  spectrum  to 
be  referred  to  simply  as  a  shock  response  spectrum. 

A '  SRS  is  a  convenient  means  of  illustrating  the  peak  response  of  a  collection  of  SDoF 
systems  to  a  particular  shock  excitation.  Due  to  its  simplicity,  the  SRS  has  become  widely 
employed  as  a  means  of  describing  the  shock  response  of  structures  and  equipment.  Years  of 
earthquake  research  have  shown  that  all  earthquake  shock  response  spectra  display  similar 
characteristics  In  fact,  it  has  been  shown  that  approximate  upper  bound  response  spectra  may  be 
constructed,  based  only  on  the  peak  displacement,  velocity  and  acceleration  of  the  oscillator  base. 
For  SDoF  systems  damped  at  between  5  and  10  percent  of  critical  damping,  an  approximate  SRS 
can  be  constructed  by  simply  multiplying  the  peak  displacement,  velocity  and  acceleration  of  the 
base  by  factors  of  1.0,  1.5,  and  2.0,  respectively.  Approximate  shock  response  spectra  generated 
by  this  approach  are  assumed  to  represent  the  upper  bound  of  the  actual  shock  response  spectra, 
and  are  further  considered  to  be  independent  of  the  precise  form  of  the  input  motion 

A  technique  similar  to  that  described  above  for  bounding  earthquake  shock  response  spectra 
has  been  adopted  by  conventional  weapon  effects  analysts.  Kiger,  et.  al.  (1984),  have  shown  that 
in-structure  shock  response  spectra  can  be  bounded  by  multiplying  the  peak  in-structure 
displacement,  velocity,  and  acceleration  by  factors  of  12,  1.5,  and  2  0,  respectively.  Shock 
response  spectra  generated  by  this  techni(]ue  are  assumed  to  give  an  upper  bound  on  the  response 
of  a  damped  (5  to  10  percent  of  critical)  oscillator  located  near  the  center  of  a  buried  facility. 

Since  in-structure  motions  are  t)pically  described  using  shock  response  spectra  (actually 
tripartite  shock  response  spectra),  it  was  only  ni  tural  to  attempt  to  quantify  equipment  failure  in 
terms  of  shock  spectra.  The  advantage  of  expressing  equipment  failure  with  the  same 
representation  used  to  quantify  the  in-structure  shock  environment  is  obvious.  Expressing 
equipment  fragility  in  terms  of  shock  response  spectra  greatly  simplifies  the  design  of  equipment 
shock  isolation 

The  approach  employed  to  determine  equipment  fragility  is  to  subject  an  item  of  equipment  to 
a  base  motion  of  a  certain  intensity,  frequency,  etc  If  the  equipment  does  not  fail,  the  intensity  of 
the  base  motion  is  increased  until  failure  occurs.  The  shock  response  spectrum  of  the  base  motion 
which  just  causes  the  equipment  to  fail  is  called  the  shock  tolerance  spectrum,  or  the  fragility 
spectrum  of  the  equipment,  even  though  the  test  SRS  is  input  waveform-dependent. 
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Figure  2.  Tripartite  Response  Spectra  for  Damping  Ratios  of  0,  2,  5  and  10  Percent 
of  Critical  (Paz,1985). 


By  comparing  the  fragility  spectrum  of  an  item  of  equipment  to  the  in-structure  SRS  at  the 
proposed  equipment  location,  shock  isolation  requirements  are  easily  assessed.  If  the  equipment 
fragility  spectrum  exceeds  the  in-structure  SRS  at  all  frequencies,  no  shock  isolation  is  required. 
If  the  in-structure  SRS  exceeds  the  equipment  fragility  spectrum  at  any  frequency,  shock  isolation 
is  required. 

This  approach  to  shock  isolation  assessment  is  oversimplified,  and  ignores  several  sources  of 
error.  The  most  important  source  of  error  is  the  fact  that  an  SRS  does  not  correspond  to  a  unique 
input  time-history.  An  infinite  number  of  different  base  motions  can  generate  a  given  SRS.  These 
different  base  motions  could  vary  greatly  in  duration,  frequency  content,  and  amplitude.  The  main 
assumption  behind  the  SRS  approach  to  equipment  fragility  is  that  equipment  failure  is 
independent  of  the  input  waveform.  In  reality,  equipment  fragility  spectra  are  probably  only  valid 
for  waveforms  resembling  the  one  under  which  the  item  of  equipment  was  actually  tested. 
Extrapolating  equipment  fragility  based  on  nuclear  testing  to  the  conventional  weapon 
environment  is  now  done  routinely  The  validity  of  this  extrapolation  has  never  been  verified 

C.  SCOPE 

This  report  quantifies  some  of  the  errors  introduced  by  assessing  the  shock  tolerance  of 
equipment  via  the  classic,  linear,  SDoF  shock  response  spectrum.  Two  simple  mathematical 
equipment  models  were  developed  for  the  investigation:  (1)  a  spring  pendulum,  which  possesses 
two  degrees  of  freedom  and  geometric  nonlinearity,  and  (2)  a  clamped  beam  carry  ing  a  tip  mass 
having  both  translational  and  rotary  inertia,  which  possesses  an  infinite  number  of  linear, 
orthogonal  modes.  These  models  were  used  to  investigate  the  effect  of  factors  such  as; 
waveform  dependency  of  equipment  shock  response,  multi-degree-of-freedom  (MDoF)  vs  SDoF 
equipment  response,  the  effect  of  multidirectional  loading  on  equipment  response,  and  the  effect 
of  geometric  nonlinearity  on  equipment  response. 
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SECTION  n 


LITERATURE  REVIEW 


A.  INTRODUCTION 

llie  literature  on  equipment  shock  tolerance  is  extensive.  It  includes  literature  on  methods  of 
analyzing  the  transient  response  of  mechanical  systems  to  shock  input,  the  origin  of  and 
assumptions  underlying  the  shock  response  sp}ectrum,  observed  equipment  damage  due  to  shock, 
equipment  shock  testing  methods,  and  equipment  shock  test  data  analysis  and  interpretation. 
Since  the  purp>ose  of  this  effort  is  to  examine  the  accuracy  of  the  shock  response  spectrum  in 
predicting  the  shock  tolerance  of  a  representative  equipment  subsystem,  an  examination  of  the 
history  of  the  shock  response  spectrum  was  in  order.  This  inquiry  led  to  a  cursory  but  fascinating 
review  of  the  theory  of  the  transient  response  of  linear  systems,  as  developed  by  well-known 
mathematicians  and  engineers  including  Cauchy,  Heaviside,  Steinmetz,  Carson,  Bush,  Biot, 
Kousner,  and  Gardner  and  Barnes. 

Heaviside  has  been  credited  with  the  development  of  much  of  modem  operational  system 
theory,  and  Bush's  book  reinforces  that  notion.  However,  as  Bush  and  colleagues,  Gardner  and 
Barnes,  studied  the  subject  further,  it  became  clear  that  Heaviside  had  ignored  more  fundamental 
methods  developed  earlier  by  Cauchy  and  other  mathematicians.  Steinmetz,  often  credited  with 
introducing  complex  numbers  into  AC  circuit  analysis,  gets  no  credit  from  Bush  or  Gardner  and 
Barnes.  Biot's  Ph  D  thesis,  written  at  Cal  Tech  in  1932  under  Von  Karman,  applied  complex 
variable  methods  to  transient,  linear  structural  dynamics,  but  the  shock  response  spectrum 
concept  is  not  clearly  evident  there.  It  is  clearly  evident  in  Housner's  1941  Cal  Tech  Ph  D.  thesis. 
It  was  Jacobsen,  Crede,  and  particularly  Newmark  who  popularized  the  shock  response  spectrum 
concept  in  protective  constuction  shock  isolation  analysis. 

B  EQUIPMENT  FRAGILITY  BASED  ON  THE  SRS  METHOD 

Current  design  guidelines  for  isolating  equipment  from  the  effects  of  conventional  weapon- 
induced  structural  motions  are  based  on  SAFEGUARD  test  data  The  SAFEGUARD  Hardness 
Assurance  Program  was  conducted  to  investigate  the  reliability  of  equipment  installed  in  hardened 
Ballistic  Missile  Defense  (BMD)  facilities  In  SAFEGUARD,  some  300  items  of  off-the-shelf 
commercial  equipment,  assigned  to  32  generic  equipment  groups,  were  tested  to  qualify  30,000 
critical  items  located  in  SAFEGUARD  installations.  Both  electical  and  mechanical  components 
were  subjected  to  shake-table  testing  The  base  motion  consisted  of  a  sine  sweep  pulse  of  5 
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second  duration,  selected  to  fit  a  prescribed  acceleration  spectrum  The  prescribed  acceleration 
spectrum  was  based  on  calculations  of  the  shock  environment  within  the  SAFEGUARD  BMD 
facility. 

Initially,  fragility  testing  was  conducted  on  selected  equipment  items,  but  most  of  the  testing 
was  limited  to  proof-type  qualification  tests,  i  e.  tests  performed  at  shock  input  levels  at  or  below 
design  levels,  to  qualify  the  equipment  for  operational  service.  Most  SAFEGUARD  data  is  given 
in  terms  of  the  base  motion  SRS  that  the  equipment  survived  in  proof-type  qualification  tests.  In 
most  cases  there  is  no  indication  of  how  close  the  SRS  is  to  the  failure  threshold  of  the 
equipment.  In  addition,  much  of  the  data  does  not  encompass  the  sensitive  frequencies  for  the 
items  of  equipment. 

Despite  the  obvious  limitations  of  the  SAFEGUARD  data,  all  land-based  equipment  shock 
isolation  designs  are  currently  based  on  this  information.  The  SAFEGUARD  procedure  is 
recommended  by  both  the  Army's  Fundamentals  of  Protective  Design  for  Conventional  Weapons 
and  the  Air  Force's  Protective  Construction  Design  Manual,  because  a  better  procedure  is  not  yet 
available. 
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SECTION  m 


EQUIPMENT  SHOCK  RESPONSE 


A.  INTRODUCTION 

The  previous  sections  have  raised  several  questions  regarding  the  validity  of  defining 
equipment  shock  tolerance  via  the  tripartite  SRS.  The  objective  of  this  section  is  to  quantify  the 
magnitude  of  the  errors  intioduced  by  the  SRS  approach  .  Two  analytical  models  were  developed 
for  this  purpose:  (1)  a  spring  pendulum,  and  (2)  a  clamped  beam  carrying  a  tip  mass. 

B,  FRAGILITY  SPECTRA  FROM  AN  IDEAL  EXPERIMENT 

The  greatest  error  associated  with  the  SRS  approach  to  equipment  shock  tolerance 
characterization  lies  in  the  assumption  that  the  shock  tolerance  of  the  equipment  is  independent  of 
the  precise  input  waveform.  To  examine  the  validity  of  this  assumption,  the  authors  designed  a 
perfect  analytical  experiment  using  an  ideal  item  of  equipment  with  known  fragility.  It  was 
assumed  that  the  item  of  equipment  is  perfectly  modeled  by  an  SDoF  undamped  oscillator  with  a 
natural  frequency  of  25  Hz.  This  simple  item  of  equipment  was  assumed  to  fail  when  the  relative 
displacement  between  the  mass  and  the  support  reached  a  value  of  0.25  inches 

Fragility  testing  of  this  simple  SDoF  equipment  item  was  conducted  analytically  The  support 
excitation  was  a  half-cycle  sine  velocity  pulse  of  the  form 

>'(0 ^ for  O^t^t^  (4) 

and 

y(/)  =  0  otherwise  (5) 

where  is  the  amplitude  of  the  pulse,  and  tj  is  the  duration. 

The  first  fragility  test  was  conducted  for  a  half-cycle  sine  velocity  pulse  with  a  duration  of  1 
second.  The  amplitude  of  the  support  motion  was  increased  until  the  maximum  relative 
displacement  reached  a  value  of  0  25  inches.  Equipment  failure  occured  when  the  amplitude  of 
the  sine  pulse  reached  a  value  of  982  in/sec. 

A  second  fragility  test  was  performed  with  a  pulse  duration  of  0.01  seconds  The  pulse 
amplitude  producing  failure  was  42  in/sec. 
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Thus,  we  now  have  two  half-cycle  sine  velocity  pulses  which  will  just  precipitate  failure  in  the 
ideal  SDoF  25  Hz  item  of  equipment.  The  maximum  relative  displacement  spectra  for  these  two 
sine  pulses  are  shown  in  Figure  3.  Notice  that  the  two  spectra  coincide  at  only  a  single  point 
(Frequency  =  25  Hz). 

The  tripartite  shock  response  spectra  for  the  two  half-cycle  sine  base  excitations  are  shown  in 
Figure  4  Again,  notice  that  the  spectra  coincide  at  only  one  point  (the  point  at  which  they  were 
forced  to  agree).  Which  of  these  two  spectra  is  the  correct  fragility  spectrum  for  the  item  of 
equipment?  Neither  one!  There  are  an  infinite  number  of  possible  fragility  spectra  for  this  simple 
item  of  equipment.  The  only  similarity  between  these  spectra  is  the  fact  that  they  all  intersect  at  a 
frequency  of  25  Hz  and  an  amplitude  of  0.25  in. 

This  analytical  experiment  clearly  illustrates  the  inadequacy  of  SRS-based  fragility 
characterization 

C  EFFECT  OF  TWO-DIMENSIONAL  BASE  MOTION  ON  EQUIPMENT  RESPONSE 

An  item  of  equipment  located  in  a  hardened  shelter  will  be  subjected  to  three-dimensional 
support  motion  In  all  likelihood,  the  equipment  item  possesses  multiple  degrees  of  freedom,  and 
exhibits  nonlinear  response  to  shock  input.  The  obvious  question  is  -  what  is  the  error 
introduced  by  modeling  a  nonlinear  MDoF  item  of  equipment  subject  to  three-dimensional 
support  motion,  with  a  linear  SDoF  oscillator? 

To  begin  a  quantitative  investigation  of  this  question,  an  equipment  model  capable  of 
responding  to  multi-directional  support  motion  was  required.  A  simple  two-dimensional  spring 
pendulum  equipment  model  was  developed  for  the  investigation  (see  Figure  5).  This  simple 
model  possesses  two  degrees  of  freedom  and  geometric  nonlinearity.  The  model  is  capable  of 
responding  to  horizontal  and  vertical  support  motion,  individually  or  in  combination.  The 
complete  mathematical  development  of  the  spring  pendulum  model  is  presented  in  Appendix  A 
The  computational  algorithm  described  in  Appendix  A  was  coded  into  a  FORTRAN  computer 
program  named  SPEND.  A  listing  of  SPEND  is  included  as  Appendix  B. 

The  spring  pendulum  equipment  model  exhibits  characteristics  of  both  a  spring-mass  oscillator 
and  a  simple  pendulum.  Thus,  it  possesses  two  characteristic  natural  frequencies  of  interest,  the 
natural  frequency  of  a  simple  pendulum,  given  by 

(H2)  (6) 

where  g  is  the  acceleration  due  to  gravity,  and  /  is  the  length  of  the  pendulum,  and  the  natural 
frequency  of  a  simple  spring-mass  oscillator,  given  by 
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Figure  4.  Shock  Response  Spectra  for  Half-Cycle  Sine  Velocity  Pulses. 


where  k  is  the  spring  stiffness,  and  m  is  the  oscillator  mass. 

A  parameter  study  was  performed  to  investigate  the  response  of  the  spring  pendulum 
subjected  to  simultaneous  horizontal  and  vertical  support  motion.  As  a  baseline  for  comparison, 
the  response  to  the  superimposed  motion  was  compared  to  the  response  due  to  vertical  base 
motion  only.  The  support  motion  for  all  cases  was  a  half-cycle  sine  velocity  pulse  of  the  form 
described  in  Equations  (4)  and  (5).  All  support  motions  had  an  amplitude  of  20  in/sec  Pulse 
durations  varied  from  0.05  to  0.5  seconds. 

The  first  set  of  calculations  were  performed  using  velocity  pulses  with  a  duration  of  0.05 
seconds.  The  results  are  given  in  Figure  6.  The  solid  line  in  Figure  6  represents  the  response  of 
the  spring  pendulum  subjected  to  vertical  support  motion  only  This  response  mode  corresponds 
to  a  simple  1-D  spring-mass  oscillator.  The  three  remaining  curves  represent  the  response  of  the 
spring  pendulum  when  subjected  to  both  horizontal  and  vertical  support  motion  simultaneously. 
Notice  that  biaxial  support  motion  always  results  in  a  spring  force  greater  than  that  caused  by 
vertical  support  motion  only.  The  amplification  of  the  spring  force  increases  as  the  pendulum 
frequency  approaches  the  frequency  of  the  support  motion. 

Figure  7  presents  the  results  of  calculation  set  number  2.  For  this  run,  the  pulse  durations 
were  all  0.1  second. 

The  results  of  the  third  calculation  set  are  given  in  Figure  8.  The  support  motions  for  these 
calculations  had  a  duration  of  0.25  seconds.  This  corresponds  to  a  frequency  of  approximately  4 
Hz  Notice  the  marked  increase  in  the  spring  pendulum  response  for  the  pendulum  frequency  of  3 
Hz. 

Results  of  the  last  calculation  set  are  given  in  Figure  9.  The  duration  of  the  support  motion  for 
this  set  was  0  5  seconds  (frequency  --  2  Hz)  Notice  the  response  of  the  spring  pendulum  for  the 
pendulum  frequency  of  2  Hz  When  the  frequency  of  the  support  mot'cn  and  the  pendulum 
frequency  of  the  spring  pendulum  model  coincide,  the  response  of  the  model  becomes  erratic. 

Figure  9  dearly  illustrates  what  is  possible  when  a  nonlinear  equipment  item  with  certain 
fundamental  frequencies  is  subjected  to  biaxial  support  motion.  The  ma'  irnum  response  of  such 
an  equipment  item  could  be  many  times  greater  than  the  response  calculated  using  a  simple 
oscillator. 
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Figure  6.  Response  of  Spring  Pendulum  to  Combined  Horizontal  and  Vertical  Support  Motion 
(Base  Input:  Amplitude  =  20  in/sec.  Duration  =  0.05  sec) 
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Figure  7.  Response  of  Spring  Pendulum  to  Combined  Horizontal  and  Vertical  Support  Motion 
(Base  Input:  Amplitude  =  20  in/sec,  Duration  =  0. 10  sec). 
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Figure  8,  Response  of  Spring  Pendulum  to  Combined  Horizontal  and  Vertical  Support  Motion 
(Base  Input;  Amplitude  =  20  in/sec.  Duration  "  0  25  sec) 


Spring-tlass  Frequency  (Hz) 

Figure  9.  Response  of  Spring  Pendulum  to  Combined  Horizontal  and  Vertical  Support  Motion 
(Base  Input:  Amplitude  =  20  in/sec.  Duration  =  0.50  sec). 


D  EFFECT  OF  MDOF  EQUIPMENT  RESPONSE 


It  is  conunon  practice  to  use  the  fundamental  period  (or  frequency)  of  an  MDoF  systen\,  and 
enter  the  shock  response  spectrum  for  a  given  support  motion  to  determine  the  peak  structural 
response  of  the  system.  The  contributions  to  the  structural  response  from  higher  modes  are 
usually  ignored. 

To  quantify  the  errors  introduced  by  ignoring  higher  mode  contributions  to  the  structural 
response,  a  cantilever  beam  carrying  a  tip  mass  having  both  translational  and  rotary  inertia  was 
used  as  a  two-dimensional  MDoF  equipment  model.  The  equipment  response  was  assumed  to  be 
linearly  elastic.  This  simple  MDoF  model  is  shown  in  Figure  10.  The  mass  density,  cross- 
sectional  area,  Young's  modulus,  moment  of  inertia,  and  length  of  the  beam  are  denoted  by  p, 

£,  /,  and  L,  respectively.  These  parameters  were  assumed  constant  throughout  the  beam  The 
mass  and  the  radius  of  gyration  of  the  tip  mass  are  denoted  by  m  and  r,  respectively 

The  derivation  of  the  equation  of  motion,  along  with  the  initiai  and  boundary  conditions  for 
this  structural  system  when  subjected  to  support  motion  are  given  in  Appendix  C  The  support 
motion  is  prescribed  as  a  displacement  time-history,  Z(t),  allowing  a  velocity  jump  at  time  t  =  0 
The  relative  displacement  of  the  beam  with  respect  to  the  support  is  denoted  as  u(x,l).  Based  on 
the  principle  of  linear  superposition,  the  total  response  of  the  system  can  be  expressed  as  the  sum 
of  modal  contributions  via  a  modal  analysis.  Due  to  the  orthogonaiity  of  different  vibration 
modes,  each  modal  equation  of  motion  can  be  solved  as  that  for  an  SDoF  system  The  derivation 
of  the  orthogonality  relationship  for  distinct  vibration  modes,  and  the  modal  equation  of  motion 
are  also  included  in  Appendix  C.  Furthermore,  the  exact  solution  can  be  obtained  if  the  support 
motion  can  be  described  by  a  simple  analytical  function. 

An  item  of  equipment  with  =  1  and  k2  =  0.0025  (refer  to  Equations  (C-51)  and  (C-52)  in 
Appendix  C),  is  used  as  an  example  The  characteristic  curve  (i  e.,  the  plot  oi  Equation  (C-58)) 
of  this  system  is  shown  in  Figure  1 1.  The  roots  of  this  curve  are  the  characteristic  values  of  the 
free-vibration  equation  (Equation  (C-27)),  and  the  corresponding  characteristic  functions 
(obtained  from  using  Equations  (C-36),  (C-49)  and  (C-59))  are  the  free-vibration  mode  shapes 
Although  this  system  ha  an  infinite  number  of  natural  vibration  modes,  only  the  first  seven  modes 
are  retained  for  a  modal  analysis  The  characteristic  values  of  these  modes  are  identified  in  Figure 
11.  The  normalized  mode  shapes  of  the  first  seven  modes  are  shown  in  Figures  12  through  18, 
respectively  The  accuracy  of  the  characteristic  values  and  mode  shapes  deteriorates  with  higher 
modes,  but  fortunately  the  contributions  from  higher  modes  are  relatively  insignificant. 

The  support  motion  (or  base  acceleration)  in  this  study  was  modeled  as  a  unit  triangular  pulse 
with  duration  (j  and  no  rise  time  Mathematically,  the  base  acceleration  can  be  expressed  as 
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Figure  12.  First  Mode  Shape  ( A,i  =  1.247 ) 


22 


.928) 


24 


Figure  15.  Fourth  Mode  Shape  (  =  8.832) 
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Normalized  position, 


igure  18. 
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where  T ,  defined  by  Equation  (C-17)  in  Appendix  C,  is  the  characteristic  period  of  the  system. 
The  effects  of  pulse  duration  on  the  equipment  response  parameters,  such  as  tip  mass 
displacement,  base  shear,  and  base  moment,  were  evaluated  as  functions  of  IT.  The  cyclic 
natural  frequency  in  Hz  of  the  ith  mode,  /j- ,  is  related  to  the  characteristic  period  by 


fi 


hL 

InT 


(9) 


The  characteristic  period,  T,  of  the  system  was  assumed  to  be  0.02  sec  in  this  example.  With 
the  characteristic  values  and  mode  shapes  given,  the  modal  amplitudes,  6,,  for  the  first  seven 


modes  were  obtained  by  directly  integrating  Equation  (C*123).  Figure  19  shows  the  relative 
displacement  time-history  of  the  tip  mass  for  a  base  motion  of  duration  2T,  considering  MDoF 
response.  The  base  shear  and  base  moment  can  be  calculated  using  Equations  (C-3)  and  (C-4), 
respectively.  Figure  20  shows  the  base  shear  time-history  for  a  base  motion  of  duration  ST. 

Since  the  response  in  each  mode  is  given  explicitly,  the  error  of  ignoring  higher  mode 
contributions  can  be  computed.  The  error  is  defined  herein  as  the  difference  between 
MDoF(including  all  seven  modes)  and  SDoF(including  the  first  mode  only)  values,  assuming  the 
total  response  of  the  equipment  is  predominantly  in  the  first  mode.  The  percentage  error  from 
using  the  SDoF  model  in  the  calculation  of  tip  displacement,  base  shear,  and  base  moment  is 
shown  as  a  function  of  impulse  duration  in  Figure  21.  The  error  increases  with  the  degree  of 
differentiation  of  the  relative  displacement  with  respect  to  length. 
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Relatiue  Displacement  of  Tip  Mass  due  to  Unit  Triangular  Impulse  at  Base 

Impulse  Duration  t^/T  =  2 
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Figure  19.  Relative  Tip  Mass  Displacement  due  to  a  Unit  Triangular 


Base  Shear  due  to  Unit  Triangular  Base  Acceleration 

Impulse  Duration  =  5 
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Figure  20.  Base  Shear  due  to  a  Unit  Triangular  Pulse  at  Base  (tjj/T-5) 
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Figure  2 1 .  Percentage  Enor  from  Using  SDoF  Model  for  MDoF  Response 


SECTION  IV 


CONCLUSIONS  AND  RECOMMENDATIONS 


A.  CONCLUSIONS 

The  objective  of  this  study  was  to  investigate  the  errors  associated  with  characterizing  the 
shock  tolerance  of  equipment  via  the  classic,  linear,  SDoF  shock  response  spectrum.  Two  simple 
mathematical  equipment  models  were  developed  for  the  investigation;  (1)  a  spring  pendulum, 
which  possesses  two  degrees  of  freedom  and  geometric  nonlinearity,  and  (2)  a  clamped  beam 
carrying  a  tip  mass  having  both  translational  and  rotary  inertia,  which  possesses  an  infinite  number 
of  linear,  orthogonal  modes.  These  models  were  used  to  investigate  and  quantify  the  effect  of 
factors  such  as.  waveform  dependency  of  equipment  shock  response,  multi-degree-of-freedom 
(MDoF)  vs.  SDoF  equipment  response,  the  effect  of  multidirectional  loading  on  equipment 
response,  and  the  effect  of  geometric  nonlinearity  on  equipment  response. 

This  study  uncovered  and  quantified  several  shortcomings  of  the  SRS  approach  to 
characterizing  equipment  shock  tolerance. 

•  It  was  proven  that  SRS-based  fragility  spectra  are  not  unique  Even  for  an  equipment  item 
which  can  be  modeled  by  a  simple  SDoF  undamped  oscillator,  different  base  excitations 
generally  produce  different  fragility  spectra. 

•  The  spring  pendulum  model  was  used  to  illustrate  the  response  of  a  geometrically  nonlinear 
item  of  equipment  to  two-dimensional  base  otion.  It  was  shown  that  extremely  erratic 
behavior  is  possible  if  the  frequency  of  the  base  motion  approaches  one  of  the  natural 
frequencies  of  the  equipment  model.  The  maximum  response  can  be  'everal  times  greater 
than  that  of  a  simple  SDoF  oscillator 

•  The  cantilever  beam  model  was  used  to  show  that  the  maximum  repsonse  of  an  MDoF  system 
will  always  exceed  the  response  of  a  simple  oscillator. 

In  conclusion,  this  study  has  proven  the  inadequacy  of  the  linear,  SDoF  shock  response 
spectrum  for  characterizing  the  shock  tolerance  of  equipment.  Complex,  nonlinear  mechanical 
equipment  subjected  to  muUidiiectional  support  motion  is  often  not  adequately  represented  by  an 
SDoF  model. 
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B  RECOMMENDATIONS 


Based  on  the  results  of  this  study,  a  more  rigorous  approach  for  assessing  equipment  shock 
tolerance  is  required  The  following  points  should  be  kept  in  mir  !  while  developing  this  new 
procedure: 

•  When  an  equipment  item  is  tested  to  determine  its  shock  tolerance,  the  test  input  waveform 
must  be  representative  of  the  anticipated  threat  input  waveform.  Multidirectional  support 
motion  must  be  reproduced.  Equipment  tests  should  excite  the  same  response  modes,  and 
produce  the  same  failure  modes,  as  the  actual  in-service  base  motion. 

•  Analytical  equipment  models  must  be  detailed  enough  to  reproduce  the  salient  features  of  the 
actual  equipment  response.  The  model  need  not  encompass  the  entire  item  of  equipment,  but 
it  must  adequately  represent  the  critical  components.  The  input  motion  to  the  model  must 
mimic  the  in-service  motion. 
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APPENDIX  A 


DYNAMIC  RESPONSE  OF  A  SPRING  PENDULUM  TO  SUPPORT  MOTION 


Dynamic  Response  of  a  Spring  Pendulum  to  Support  Motion 


Consider  the  response  of  the  spring  pendulum  shown  below,  to  support  motion. 


Y 

The  position  vector  of  the  support  point,  S,  is 

r  =  ucj  +  ve2 

and  the  position  vector  of  the  suspended  mass,  m,  is 

^  =  x?,+(/o  +  >')ej 

where 

4  =  unstretched  spring  length 
u  =  X-displacement  of  the  support 
V  =  Y-displacement  of  the  support 
X  =  X-displacement  of  m  from  the  undeformed  position 
y  =  Y-displacement  of  m  from  the  undeformed  position 

The  spring  vector  is 

K  -  f  =  (X - u)e, -t- [/q v)]e2 

so  that  the  spring  length  is 

and  a  unit  vector  pointing  along  the  spring  from  support  to  mass  is 


^  {x-u)e^+[lo-i-{y-v)]e2 


1  he  spring  elongation  is 


A  =  /-/o 


and  if  the  spring  stiffness  is  k,  the  spring  restoring  force  is 


(A-1) 

(A-2) 


(A-3) 


(A-4) 


(A-5) 

(A-6) 
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The  gravitational  force  is 

P,  =fnge^ 

The  equation  of  motion  for  the  mass  is 

/)  +  Fj  =  mh 


or 


The  static  solution  is 


ysr  ~  ^ST  ~^o 


rug 

k 


The  static  solution  defines  the  initial  conditions  for  a  dynamic  problem.  To  initiate  a  finite 
difference  solution,  we  assume  that  x  and  y  are  constant  during  the  first  time  interval. 


Therefore 


= 


(A-7) 

(A-8) 

(A-9) 

(A-10) 

(A-11) 

(A.12) 

(A-13) 

(A-14) 

(A- 15) 


(A- 16) 
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(A-17) 


>1  =ysr  + 


yoW 


Equation  (A- 10)  thus  takes  the  form 


£ 

m 


so  that 


Isrl 


Xo(A/)' 


-I/, 


(A-18) 


IHi 

f  =  =  *^1 

6/ir  ^  6 

and  Equation  (A>I  I)  takes  the  form 


so  that 


>'o  =  S 


m 


y’sT 


6 


Vo  = _ m - =: - ^^1 - 

1+k^  ^ 

6m  k  6 

From  this  point  on,  the  finite  difference  solution  proceeds  smoothly. 
Step  1 :  Xq  =  0 


Step  2; 


>'o 


f»g 

k 


Step  3 ; 


= 


g(AO^ 

61 


ST 


1  + 


gW 


61. 


ST 


(A-19) 


(A-20) 


(A-21) 
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Step  4; 


Vl+^o 


For  n ^  1 

Step  5: 

Step  6: 

Step  7 : 

Step  8: 

yn^g-M-^(yn-^'n)] 

Step  9; 

x„.,=2x„-x„_^+x„(Aty 

Step  10: 

Return  to  step  5 

APPENDIX  B 

SPEND  PROGRAM  LISTING 
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PROGRAM  SPEND 
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ooo  n  oooon  o  onn  ooo  nno  ooo  non  non 


C  *  * 

^  Hi,******************************************************** 

C 

READ(15, *)  PWT, STIFF, ULENGTH,NPTS,DELT 
PRINT  *, '  Reading  Horizontal  Base  Motion' 

READdS,*)  (U(N),  N=0,NPTS-1) 

PRINT  *, '  Reading  Vertical  Base  Motion' 

READ(15,*)  {V(N),  N=0,NPTS-1) 

C 

PRINT  ’  Computing  Pendulum  Response' 

ACCELERATION  OF  GRAVITY  {IN/SEC**2) 

AGRAV=386.0 

CALCULATE  THE  MASS  OF  THE  PENDULUM  (LB-SEC**2/IN) 

PMASS=PWT/AGRAV 

CALCULATE  PENDULUM  DISPLACEMENTS  0  TIME  0  (I.E.,  STATIC  DEFLECTION) 

X(0)=0.0 
Y{0)«=PWT/STIFF 
STATLEN=Y (0) +ULENGTH 

CALCULATE  DISPLACEMENT  0  TIME  STEF  1 

CONl= ( AGRAV*  DELT*  *  2 . ) / ( 6 . *  STATLEN ) 

X{l)=(CONl/(l+CONl) )*U(1) 

CON2= (STIFF* DELT** 2. )/ (6.*PMASS) 

Y(1)=Y(0)+(CON2/(1+CON2) )*V(1) 

INITIALIZE  ACCELERATIONS  &  SPRING  FORCE  0  TIME  0 

SF{0)=PWT 
XDD(O) =0. 0 
YDD(O) =0. 0 

LOOP  THRU  TIME  STEPS  FROM  N=1  TO  NPTS-1 
DO  100  N-1, NPTS-1 

SPLEN=DSQRT( (DABS (X(N) -U(N) ) ) * *2 . + (ULENGTH+ ( Y (N ) -V (N) ) ) * *2 . ) 

ERROR  TRAP:  FOR  LOW  FREQUENCY  SYSTEMS,  THE  SPRING  MAY  BECOME 
FULLY  COMPRESSED.  THIS  LEADS  TO  A  DIVIDE  BY  ZERO 
ERROR  WHEN  CALCULATING  BETA. 

IF  (SPLEN  .LT.  (0. 1*ULENGTH) )  THEN 

PRINT  *, '  Excessive  spring  compression  has  occured!' 

STOP 
ENDIF 

BETA= (STIFF/PMASS) * ( 1-ULENGTH/SPLEN) 

CALCULATE  ACCELERATIONS  0  TIME  STEP  N 

XDD(N)=-BETA* (X(N)-U(N) ) 

YDD(N)=AGRAV-BETA* (ULENGTH+ ( Y (N ) -V (N) ) ) 

C 


■ 
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onn  nan  n  n  n  non  ooooooo  OOOOO 


SPRING  FORCE  0  STEP  N 


SF (N)  =  ( SPLEN-ULENGTH ) *  STI FF 

DISPLACEMENT  0  STEP  N+1 

X{N+1) =2 . *X (N) -X(N-l) +XDD(N) *DELT**2 
Y(N+1)-2.*Y(N)-V(N-1) +YDD(N) *DELT**2 


100  CONTINUE 
WRITE  OUTPUT  FILES 


PRINT  Saving  Penduluin  Response 

X-DISPIACEMENT 


WRITE (16,*)  NPTS 
WRITE (16,*)  DELT 
WRITE(16,*)  {X(N), 

CLOSE(16) 

Y-DISPLACEMENT 

WRITE(17,*)  NPTS 
WRITE (17,*)  DELT 
WRITE(17,*)  (y(N), 
CLOSE (17) 

SPRING  FORCE 


N=0,NPTS-1) 


N=0,NFTS-1) 


WRITE(18,*) 
WRITEde,  *) 
WRITEde,  *) 
CLOSE (18) 

NPTS 
DELT 
(SF(N) , 

N=0,NPTS-1) 

X-ACCELERATION 

WRITE(19, «) 
WRITEd9,  ♦) 
WRITEd9,  •) 
CLOSE (19) 

NPTS 

DELT 

(XDD{N) 

,  N=0,NPTS-1) 

Y-ACCELERATION 

WRITE (20, *) 
WRITE (20,*) 
WRITE(20,  *) 
CLOSE (20) 

NPTS 

DELT 

(YnD(N) 

,  N=0,NPTS-1) 

STOP 

END 
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APPENDIX  C 

TRANSIENT  RESPONSE  OF  A  CANTILEVER  BEAM  CARRYING  A  TIP  MASS 
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Tnnslent  Response  of  a  CantUever  Beam,  Carrying  a  Tip  Mass  Having  Both  Translational 
and  Rotary  Inertia,  to  Support  Motion 


Consider  the  cantilever  beam  shown  below,  carrying  a  tip  mass  having  both  translational  and 
rotary  inertia. 


V 


1.  Governing  Equation.  Initial  and  Boundary  Conditions 
The  total  beam  deflection  is 

y-z  +  u  (C-1) 

where 

2  =  support  motion 
u  -  beam  deflection  WRT  the  support 
The  equation  of  motion  for  a  beam  element  is 

- dx-(Adx — ^^qAdx — r(z  +  u)  (C  2) 

dx  dr  dr 

Neglecting  the  beam  element  rotary  intrtia  yields 
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and  the  beam  moment-curvature  relation  is 


M=n- 


Substituting  Equations  (C-3)  and  (C-4)  into  Equation  (C-2)  yields 


(C-3) 


(C-4) 


d*u  pA  d^u  _  pA  d^z 


The  initial  and  boundary  conditions  are: 


Now  let 


so  that 


(C-5) 


>(r,0)  =  0 

(O^x^L) 

(C-6) 

|(..0)  =  0 

(O^r^I) 

iC-1) 

«(o,o  =  o 

(/^O) 

(C-8) 

1^(0./)  =  0 
dr 

(r^O) 

(C-9) 

£/0(i,/) 

II 

1 

Vis'' 

(/SO) 

(C-10) 

"'0 

(/SO) 

(C-ll) 

u 

a  =  — 

L 

(C-12) 

H 

(C-13) 

(C-14) 

d  _ 

_  d  W 

(C-15) 

dr 

''  dl,dx~  Ld^ 
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Then  Equation  (C-5)  can  be  written  in  the  form 


d*a  pAL*  d^a  _  pAL* 
d^*  Ef  d/^  '  Ef  dt^ 

and  if  we  set 


/ 


so  that 


^  _  9  (Ef/  _  I  d 
9/  d\f  dt  7"  9v 

then  Equation  (C-16)  can  be  written  in  the  final  dimensionless  form 


d*a  d^a  _  9^p 
9V  d\if^  ^ 


with  the  initial  and  boundary  conditions 

aa,0)  =  0  (O^^Sl) 


|^«.0)  =  -|£(0) 

9\jr  9v 


(C-16) 

(C-17) 

(C-18) 

(C-19) 

(C-20) 

(C-21) 

(C-22) 


a(0,v)  =  0 
~(0,V)  =  0 


,  mr 


(V^O) 

(vso) 


—  (l.vr)  =  -- 
9^'  pAL^ 


m 


37r(l.  V')  =  “77 
95  pAL 


d^a 


ilW) 


We  first  consider  the  homogeneous  form  of  Equation  (C-20), 

9^a  9^a  „ 


(V^O) 

(V*/S0) 


(C-23) 

(C-24) 

(C-25) 

(C-26) 

(C-27) 
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because  it  turns  out  to  generate  a  family  of  orthogonal  functions  (free  vibration  modes)  that  can 
be  used  to  solve  the  nonhomogeneous  problem 
We  assume  that 


(C.28) 


so  that,  setting 


de 

d\(f 


=e 


substituting  Equation  (C-28)  into  Equation  (C-27)  yields 

0+^  =  0 

or,  dividing  by  ({>6, 


(C-29) 

(C-30) 

(C-31) 


e 

■2— +- 

<i>  e 


=  0 


(C-32) 


Equation  (C-32)  says  the  sum  of  two  independent  quantities  (one  a  function  of  space,  the  other  a 
function  of  time)  is  zero;  therefore  they  must  both  be  constant.  We  therefore  set 


(p  e 

so  that  ^  and  6  must  obey  the  equations 

-X‘'<|>  =  0 

0+ve=o 

The  solution  to  Equation  (C-34)  can  be  written  in  the  form 

<(.=  /!/,+ 5/, +C/3  +  D/, 

where 

/,  =  cosh  +  cosX^ 

/j  =  sinh  X4  -  sin 
/,  =  cosh  cos 
/,  =  sinh  +  sin 

so  that 

<l.'  =  X(/l/,  +  5/,  +  C/,+D/,) 
<f"  =  X^(>f/3  +  5/,  +  C/.+D/,) 
4>"'  =  ?^’(-4/,  +  5/,+C/3+D/,) 


(C-33) 

(C-34) 

(C-35) 

(C-36) 

(C-37) 

(C-38) 

(C-39) 

(C-40) 

(C-41) 

(C-42) 

(C-43) 


55 


V  =  ^(4/-,  +  B/.  +C/,  +D/J  =  V* 

Now 

/,(0)  =  2 

/:(0)  =  /3(0)  =  /,(0)  =  0 

so  that  Equations  (C-23)  and  (C'36)  yield 

<KO)  =  2y4=0 

and  Equations  (C»24)  and  (C-41)  yield 

<t»'(0)  =  2U)  =  0 

so  that 

A  =  D^Q 

If  we  now  set 

/.0)  =  f\ 


pAL 


(iJ- 


then  Equations  (C-25),  (C*41),  and  (C-42)  yield 

or 

X'  {BF,  +  CF, )  =  k,k^\\BF,  +  CF, ) 

or 

5(F,  -Jt,JtjX’F3)+C(f;  -it.itjVF,)  =  0 
and  Equations  (C-26),  (C-36),  and  (C-43)  yield 

r{\)e=k,me 

or 

V  (5F,  +  CFj )  =  -k,\*{BF,  +  CF, ) 

or 

5(F;  +  k.XF^ )  C(Fj  +  )  =  0 

Equations  (C*54)  and  (C-56)  yield  the  relative  values  of  B  and  C. 


C  ^  F,  -  ^  +  /t.XF;  ^ 

F  "  -  Jt./tjX’F,  "  Fj  +  k,KF,  ■  ^ 

so  that 

{F,^k}S,){F,-k,k,\^F,)-U'\+k}J^,){F,-k,k;K^F,)  =  Q 


(C-44) 

(C-45) 

(C-46) 

(C-47) 

(C-48) 

(C-49) 

(C-50) 

(C-51) 

(C-52) 

(C-53) 

(C-54) 

(C-55) 

(C-56) 

(C-57) 

(C-58) 
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For  each  root  of  Equation  (C-58),  X.,,  there  is  a  ratio 


given  by  Equation  (C*57). 


(C-59) 


2.  Orthogonality  of  Vibration  Mode  Shapes 

Now  consider  the  incremental  work  done  by  the  end  actions  associated  with  mode  i,  acting 
through  incremental  displacements  associated  with  mode  j. 


(C-60) 

Now 

ly  r-r  ^  ^  ^ 

(C-61) 

ax’ 

(C-62) 

c 

II 

(C-63) 

ax 

(C-64) 

so  that  Equation  (C-60)  can  be  written  in  the  form 

5»'.  =  fe.5e,[C4';-0,I 

(C-65) 

L 


y  o,5e,0o: -o:k 


Eliminating  between  Equations  (C-65)  and  (C-66)  yi<;lds 


(C-66) 
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Since  the  LHS  of  Equation  (C*67)  is  symmetric  with  respect  to  i  and  j,  the  RHS  must  be  also,  so 
that 


(C-68) 

Eliminating  between  Equations  (C-67)  and  (C-68)  yields 

(a‘,  -  A* *[*x-  m  -  o. I = 0 

(C-69) 

Now 

-e- 

o 

II 

o 

(C-70) 

<})'(0)  =  0 

(C-7I) 

and  if  we  set 

(}K1)  =  <1) 

(C.72) 

then  Equations  (C-33)  and  (C-53)  yield 

<2)"  =  X,XjX^O' 

(C-73) 

and  Equations  (C-33)  and  (C-55)  yield 

O'"  =  -/t,X*<D 

(C-74) 

so  that  Equation  (C-69)  yields 

(K  +*,A'<i.,o,] 


=  (A‘,  -  =  0  (C-7S) 

When  j  *  i  the  second  factor  on  the  LHS  of  Equation  (C-75)  must  be  zero,  and  when  j  =  i  that 
factor  is  arbitrary.  Thus  we  can  set 

=  5,,  (C-76) 

Equation  (C-76)  would  be  most  useful  if  the  integral  were  expressed  in  terms  of  boundary  values. 
This  can  be  done  by  returning  to  Equation  (C-69),  and  writing  it  in  the  form  [cf  Prescott 
(1924  223),  Eq.  (9,111)] 


[♦;>.-i'>:7't'>r-»,4'rll 

Treating  X,  as  a  constant  and  as  a  variable,  applying  L'Hospital's  rule  to  the  case  in 
which  X^  =  X,,  and  noting  that 


(C-ll) 
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so  that 


Equation  (C-77)  yields 


<!>'  = — f<|>(^)l  =  — =  X 


{^[<'r'^.-o:+(0:)^-4';c]}' 


4a: 

4 

Therefore,  when  /  =  y,  Equation  (C-76)  can  be  written  in  the  form 

<I>,'+2;t,cl),<D;  + it  (<!>:)'  ,  ,  ,, 


(l  +  4A,  )<D,^  +  2A,c1),<d;  +  A,Jt,(4  +  A,A,aO(o:)'  =  4 

Now  Equations  (C-  36),  (C-37),  a’  -  t,C-50)  yield 

<P  =  BF^+Cf\ 
cD'  =  X(^/'',  +  C/';) 

and  therefore 

<D'  =  +  2F,F,BC  +  F^' 

<ixX>'  =  X\fF,B'  +{F,F,  +  F/)5C  +  F,F,C] 

(o')'  =  A^(/*/^*  +  2F,F,BC  +  F/C') 
and  thus  Equation  (C-82)  can  be  written  in  the  form 

A,F-  A,FC+A,C’  =4 

where 


(C-78) 

(C-79) 

(C-80) 

(C-81) 

(C.82) 

(C-83) 

(C-84) 

(C-85j 

(C-86) 

(C-87) 

(C-88) 


A,  =  {l  +  Ak,)F,^+2kXF,F,  -^k.k.iA  +  k.kXWA^ 

Aj  =  2[(l  +  4A:,)FjF3  +Ar,A,(FjF,  +  Fj*)  +  *,*3(4  +  A,AjA^)A’FjF,] 

Aj  -  (1  +  4it, )F/ +  2A,A,^j^4  +  W 

Substituting  Equation  (C-59)  into  Equation  (C-88)  then  yields 

b’(a,+a,p,+a,p;)=4 

$0  that 


(C-89) 

(C-90) 

(C-91) 


(C-92) 


B  = 

V^i  +  A^p.  +  Ajp? 


(C-93) 


C,=PA 

Finally,  combining  Equations  (C-67),  (C-73),  (C-74),  and  (C-76)  yields 


(C-94) 


(C-95) 


m  of  Motion 


?n  s  Pnncii 


v(x,i)  =  ^r\,ix)q,it) 


K0  =  «,(0 


>'(x,/)  =  w(r./)  +  z(/) 


=  Xn,(»H(0+w,(0 


The  kinetic  energy  of  the  system  is 


and  the  strain  energy  of  the  system  is 


2  •'«'  V  ^  ; 


(C-96) 


(C-97) 


(C-98) 


60 


Applying  Hamilton's  Principle  to  obtain  the  equation  of  motion: 


where 


<^>1 

8T  =  j^AySydx  +  my^Sy^  +mr^u[Su[ 
SU  =  {nu^Su^dx 

Jo 


^  =  &j  =  '^Tj,Sq, 
.Vi 

>■1 


<5A  =  &i =Xn.(^)5?, 

Ib'i 


!•! 


“I  =£n,'(i)4 


—I 


I«I 

Equation  (C*99)  can  be  expressed  as 

j''l^j^pA('^riq  +  uj{^r!5q}dx  +  m('^ri(L)q  +  uJ{j^Ti(L)Sg) 


(C-99) 

(C-100) 

(C-101) 

(C-102) 

(C-103) 

(C-104) 

(C-105) 

(C-106) 

(C-107) 

{C-108) 

(C-109) 

(C-110) 

(C-111) 
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■^mr 


'(X  I'CiwXl  riU)Sq)  -  j;  £/(X  n'<7XX  rt"  &l)d,y  =  0  (C-l  12) 


f  P'*(X  ')?+"J(X  n^)‘* = XX?.^/ 0’'* ’),'i,<*+",X*j.|/’'^  lA 

l.l  ;-l  <-l 

4,'^niL)q  +  Uf)('^lKL)Sq)  =  m^'^q,5q^n,(L)l)jiL)  +  mu^'^6qML)  (C-l  14) 


!•!  /•) 


ml 


I»1  J«1 


r  ^"^)(X  n"Sq)dx  =  XX^.^; 

.»1  ;«1 

Integrating  Equation  (C-l  12)  by  parts  and  using  Equations  (C-l  13)  through  (C-l  16), 

X Xf \jP^  I  - r 1 

Kl  y.lV.  I'l  ‘  / 

1=1  V  "t  ‘ 


+'”X xf  n, ( ^) H; (^)<7,&7; - J,’’ ^.Sq,r].(QnA^)^'^ 

I«1  ju\  ^ 

+'”Xf",n,(^)5<7,['  -  f 

1=1  '  '  •' 

'  X X ( a'( ( ^)9, 5?y  [|  - 1'’  5?;  n.'( ^ )  n;  ( ^ ] 

f«l  y=l  ' 


1*1  ;■! 

Rearranging  terms  yields 


(C-l  15) 


(C-l  16) 


(C-l  17) 


XXi:{4j>n,n  ^dx  +  mT],i^)nM)-i-  mr^ri'XL)T]]{L)^  +  q,  pJn',i]''dx^5q  ,dt 

1*1  j»\  ‘  ^  ^ 
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(C-118) 


=  f n,dx-k-mnXL'^Sq,dt 
The  orthogonality  relations  have  been  sho\vn  as 

j^A  T], Tj^dic +mTi,(L)Tij(L)+mr^rj;(L) rj' (L)  =  (C- 1 1 9) 

and 


Jbn,'V]>=0 


(C-120) 


for  i  *  j  j 

for  /  =  y 

Substituting  Equations  (C-1 19)  and  (C-120)  into  Equation  (C-118), 

Since  hq,  are  arbitrary  variations,  it  is  necessary  that 

for  all  f%.  Equation  (C-1 22)  is  the  equation  of  motion  for  mode  i.  The  nondimensional  form  of 
Equation  (C-1 22)  can  be  expressed  as 


(C-121) 


(C-1 22) 


e,+y9, 


(C-1 23) 


63 


